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A POSSIBLE HIDDEN POPULATION OF SPHERICAL 

PLANETARY NEBULAE 


Noam Soker and Eyal Subag 

ABSTRACT 

We argue that relative to non-spherical planetary nebulae (PNs), spherical 
PNs are about an order of magnitude less likely to be detected, at distances of 
several kiloparsecs. Noting the structure similarity of halos around non-spherical 
PNs to that of observed spherical PNs, we assume that most unobserved spherical 
PNs are also similar in structure to the spherical halos around non-spherical 
PNs. The fraction of non-spherical PNs with detected spherical halos around 
them, taken from a recent study, leads us to the claim of a large (relative to that 
of non-spherical PNs) hidden population of spherical PNs in the visible band. 
Building a toy model for the luminosity evolution of PNs, we show that the 
claimed detection fraction of spherical PNs based on halos around non-spherical 
PNs, is compatible with observational sensitivities. We use this result to update 
earlier studies on the different PN shaping routes in the binary model. We 
estimate that ~ 30% of all PNs are spherical, namely, their progenitors did not 
interact with any binary companion. This fraction is to be compared with the 
~ 3% fraction of observed spherical PNs among all observed PNs. From all PNs, 
~ 15% owe their moderate elliptical shape to the interaction of their progenitors 
with planets, while ~ 55% of all PNs owe their elliptical or bipolar shapes to the 
interaction of their progenitors with stellar companions. 

Subject headings: planetary nebulae: general — stars: AGB and post-AGB — 
stars: mass loss 


1. INTRODUCTION 

Planetary nebulae (PNs) are the descendants of intermediate mass stars (initial masses ~ 
1 — 8 M 0 ). Much of the star’s initial mass is expelled during the ~ 10^ yr that the star spends 
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on the upper asymptotic giant branch (AGB). This AGB star wind usually consists of a more 
or less spherically symmetric outflow at rates of 10 ^ — 10 ^MqJT h Most PNs, though, 
possess a global axisymmetrical rather than a spherical structure in their inner region. This is 
most evidence in the spherical halos observed around many axisymmetrical PNs (e.g., Gorradi 
et ah 2003, 2004). In axisymmetrical we include also point-symmetric PNs, which are dehned 
as PNs having more than one axisymmetric substructure, but the different symmetry axes 
intersect in a common point. This class includes also the quadrupolar PNe defined by 
Manchado et ah (1996b). These axisymmetrical PNs with spherical faint halos suggest 
that only close to the termination of the AGB, and/or later in the early post-AGB phase, 
the outflow acquires the axisymmetrical structure. The axisymmetrical inner structure of 
these PNs are much brighter than the spherical halo, implying that the mass loss rate which 
formed these parts was much higher than the regular AGB wind. As pointed out in earlier 
papers, e.g.. Sober (2000; 2002), there is a positive correlation between the onset of this 
hnal intensive wind (FIW; termed also superwind) and the transition from spherical to 
axisymmetrical mass loss. This correlation is well explained by the binary model for the 
formation of non-spherical PNs (Sober 2000; 2004), where the binary can be a stellar or a 
substellar object. 

Different classihcation schemes have been suggested over the years of the different mor¬ 
phological PN types, e.g., Manchado et ah (1996a ; 2000). We will follow Sober (1997) in 
a classification scheme which is motivated in part by theoretical considerations, and which 
contains four basic types of structure. 

1. Spherical PNs. These are dehned as PNs which shows no axisymmetrical structure at 
all. The properties of these PNs were studied by Sober (2002). One prominent property 
is that these PNs loob libe the spherical halos detected around some elliptical PNs, in 
that there is no indication for a FIW (superwind). Note that the dehnition of spherical 
PNs assumed here is quite different from that of Manchado et ah (2000). Manchado 
et ah hud that according to their dehnition ~ 25% of all observed PNs are spherical, 
while Sober (1997) hnds that only ~ 3 — 4% of all observed PNs are spherical. This 
diherence in dehnition is summarized on one side by Sober (2002), and on the other 
side by Stanghellini et al. (2002). There are 14 such PNs (those with no signature of 
a FIW) in the list of 18 PNs in Table 2 of Sober (1997). 

2. Elliptical PNs. These are PNe having a large scale elliptical shape, with a shallow 
density variation from equator to poles. They may possess jets. Examples are most 
PNs listed in Table 5 of Sober (1997). 

3. Extreme-elliptical PNs. These are dehned as PNe having a large equatorial concentra¬ 
tion of mass, e.g., a ring, but no, or only small, lobes. They may possess jets. Examples 
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are most PNs listed in Table 4 of Soker (1997). 

4. Bipolar PNs. Bipolar (also called “bilobal” and “butterfly”) PNs are defined (Schwarz 
et al. 1992) as axially symmetric PNs having two lobes with an equatorial waist 
between them. Examples are in Corradi & Schwarz (1995) and most PNs listed in 
Table 3 of Soker (1997). 

Efforts have been made over the last decade to estimate the fraction of PNs belonging to 
each group, and from that to better understand the shaping of PNs. The major question is 
whether the formation of axisymmetrical PNs requires the presence of a binary companion, 
or whether a single star can do the job. We here follow earlier works of Soker (1997, 2001a,b) 
connecting different axisymmetrical structures to different types of interaction with stellar 
and substellar companions. This is summarized in section 2, where an updated estimate 
of the fraction of PNs shaped by planets is given. Motivated by new statistical works on 
extra-solar planets, and the new work of Corradi et al. (2003) on halos around PNs, we 
set a goal to update the statistics of the different classes of PNs and their relation to the 
progenitor type. Using the results of Corradi et al. (2003) with our two basic assumptions, 
we estimate in section 3 the detection fraction of spherical PNs relative to that of non- 
spherical PNs. In section 4 we present a phenomenological toy model, that shows that the 
present detection limits of faint PNs might be compatible with the presence of a large hidden 
spherical PN population. By hidden population we refer to population of PNs which are 
below the threshold for being discovered by most visible-light surveys. We limit ourselves 
to the visible band because all large morphological samples to which shaping theories are 
compared are based on visible-light images. Some PNs in the visible-hidden population 
might be still detected in the IR, but this will not affect our claims regarding the present 
morphological classihcation to which theories are compared. In section 5 we discuss and 
summarize our main results and predictions. 


2. STATISTICS AND BINARY INTERACTION 

Only ~ 60% of all systems are binary stellar systems (Duquennoy & Mayor 1991), and 
only in a fraction of these the companion is close enough to influence the mass loss from 
the PN progenitor. The fraction of observed axisymmetrical PNs, on the other hand, is 
much larger. Soker (1997) hnds from classifying the morphology of 458 PNs this fraction to 
be 96%, namely, 4% of all observed PNs are spherical. Soker (1997) considered selection 
effects, and took the true fraction of spherical PNs to be 10%. The shaping of the non- 
spherical (axisymmetrical) PNs was attributed to three basic types of interaction as follows: 
(1) Progenitors that interact with close stellar companions that avoided a common envelope 



phase, or enter to a common envelope at a very late stage of the FIW phase, and form mainly 
bipolar PNs; ~ 11% of all PNs. (2) Progenitors that interact with stellar companions via 
a common envelope phase and form mainly extreme elliptical PNs; ~ 23% of all PNs. (3) 
Progenitors that interact with substellar, i.e., planet and brown dwarf, companions via a 
common envelope phase and form elliptical PNs (with moderate or small departure from 
sphericity); ~ 56% of all PNs. 

The interaction with the substellar companion scenario had three major predictions. 

[i) Many planets will be more massive and closer to their parent stars than Jupiter is (Sober 
1994, 1996). 

{ii) For many stars to engulf a planet at their late evolutionary stages with a high probability, 
several substellar objects must be present in most of the systems (Sober 1996). 

{Hi) About half of all progenitors of PNe have planetary systems with massive planets (Sober 
1997). This prediction was based on the notion that singly evolved stars rotate too slow when 
they become AGB stars, and hence have spherical mass loss. 

The hrst two predictions were made before extrasolar planets were found, and we adopt 
them here as well. 

In Sober (1996) massive planets were regarded as planets of mass Mp > Mj, where 
Mj is Jupiter mass, while in Sober (2001b) the minimum allowed mass to cause departure 
from spherical mass loss was reduced to Mp > O.OlMj. This reduced minimum planetary 
mass was motivated by the hnding that only ~ 5% of solar-libe stars have Jupiter libe planets 
around them (Marcy & Butler 2000), and is based on non-linear effects in the AGB envelope, 
such that slow rotation caused by low mass planets can lead to axisymmetrical mass loss. 
These effects might be excitation of non-radial pulsations, or local (but not global) dynamo 
formation of magnetic cool spots; the lower temperature above the cool spots enhanced dust 
formation, hence mass loss rate. Other non-linear processes based on locally strong magnetic 
helds are possible. For example, leabage of p-modes (sound waves) were recently proposed 
as the mechanism for the formation of chromospheric spicules in the Sun (De Pontieu et al. 
2004); spicules are dynamic jets propelled upwards from the solar surface at much below 
escape speed. In AGB stars such process might enhance dust formation, and hence increase 
mass loss rate. If the magnetic spots are concentrated in the equatorial plane (Sober 2001b), 
the enhanced mass loss rate will come with axisymmetric mass loss geometry as well. 

Sober (2001a) studied the formation of (moderate) elliptical PNs by stellar companions 
orbiting the progenitor at large orbital separations. The companion accretes mass, forms 
an accretion disb, and blows two jets. This process was found to shape ~ 5 — 20% of 
axisymmetrical PNs, hence reducing the percentage of PN shaped by planets from ~ 50% 
to ~ 30 — 45%. The estimated numbers, in the stellar and substellar shaping model, in 2001 



(Soker 2001a,b) were as follows (Table 1; for the different types of the binary interaction see 
Soker 2004): (1) AGB stars which interact with a close stellar companion that avoided a 
common envelope: ~ 15%. (2) AGB stars which interacted with a stellar companion via a 
common envelope: ~ 25% (Yungelson et ah 1993; Han et ah 1995). (3) AGB stars which 
interacted with a substellar (mainly planets) companion massive enough to cause observable 
axisymmetrical structure: ~ 35%. (This is found to be ~ 15 — 20% in the present paper). 
(4) AGB stars which interacted with a wide stellar companion: ~ 15%. (5) AGB stars which 
did not interact with any companion and formed spherical PNs: ~ 10%. (This is found to 
be ~ 25 — 30% in the present paper). 

As it turned out in recent years, the previously estimated percentage of PN progenitors 
interacting with planets is higher than the observed percentage of stars hosting planets. San¬ 
tos et ah (2004), for example, hnd that ~ 3% of stars with solar-metallicity have planets. The 
number rises to ~ 25 — 30% for stars with twice solar metallicity ([Fe/H]> 0.3). Lineweaver 
& Grether (2003) extrapolate from the detected parameters space of planetary system prop¬ 
erties to that below detection sensitivity, and conclude that the real number of planetary 
system is much larger. The fraction of solar-like stars having planets of Mp > O.lMj and 
that are close enough to interact with their parent stars, is > 20%. Over all, therefore, 
~ 25% of stars have planets that are massive and close enough to influence their mass loss 
process; a fraction higher than the one deduced by Livio & Soker (2002) who used a much 
stronger constraint on planets to influence their parent stars mass loss process. 

Two effects should be considered here. First, if most PNs are formed from more massive 
than the Sun stars, which on average are more metal rich, then the percentage of progenitors 
interacting with planet increases. Second, as note by Soker (2001b), most known Sun-like 
stars that have planets around them will not form PNs at all. Their planet(s) will deposit 
enough orbital angular momentum and energy to cause most, or even all, of the envelope 
of these stars to be lost already on the red giant branch. These stars will not evolve along 
the AGB, or evolved only along the early AGB. In any case, they will have too low mass 
envelope to form a detectable PN. This reduces the fraction of PN progenitors interacting 
with planets. 

The main conclusion form this discussion is that with present statistics of extrasolar 
planets, the percentage of PN progenitors that are expected to be influenced by planets is 
much less than ~ 35%, and it is more like ~ 15 — 20%. The percentages given above for 
the different interaction types with stellar companions are based on a more robust stellar 
population studies, and hence are more secure. In addition, recent observations suggest that 
the central star of many PNs have companions (De Marco et al. 2004; Hillwig 2004; Pollaco 
2004). The lower observe percentage of stars hosting planetary systems than the percentage 



given in previous studies motivate us to examine the possibility that the true fraction of 
spherical PNs is much larger than the one observed. 


3. SPHERICAL HALOS AND MASS LOSS HISTORY 

The FIW (superwind) at the end of the AGB is preceded by lower mass loss rate episodes, 
which in most cases known are spherical. Therefore, from theoretical considerations we 
expect all PNs to have a faint halo around them. However, the detection fraction of faint 
halos is much lower. We also expect that spherical PNs will be PNs formed by AGB stars 
which did not experienced the FIW (superwind). We therefore take the following to be the 
basic, quite reasonable, assumption (Soker 2002): Spherical PNs are similar in structure to 
the faint large halo of elliptical and extreme elliptical PNs. For example, the spherical PN 
A30 is similar to the halo of NGG 6826. 

Gorradi et ah (2003) list the 49 PNs known to have halos. We estimate the total halo 
detection rate as follows. For the sample of all well resolved PNs we use the statistics of 
Soker (1997), were 4 classes are dehned (previous sections). We do not consider the spherical 
PNs (Table 2 in Soker 1997), and the bipolar PNs having large lobes and a narrow waist, as 
in general they were not observed by Gorradi et ah (2003). The rest are 275 elliptical PNs 
and 113 extreme elliptical PNs. Around these PNs detection of halos is easier than around 
the bright and highly non-spherical bipolar PNs. Therefore, the detection fraction of halos is 
X ~ 49/(275 + 113) = 0.13; if we include all PNs studied by Soker (1997), the halo detection 
fraction is y ~ 49/458 = 0.11. The random detection fraction is even lower, because Gorradi 
et ah (2003) were aiming at specihc PNs in their search for halos. 

If we use our basic assumption in this section, that spherical PNs are similar to faint 
halos around elliptical and extreme elliptical PNs, we deduce that the detection fraction of 
spherical PNs is y ~ 0.12 times that of non-spherical PNs. This implies that the fraction of 
spherical PNs is ~ 20 — 30%, rather than the observed ~ 3%. 

We stress that the detection fraction of PNs with halos as taken from Gorradi et al. 
(2003) is an upper limit. This is because the search was aiming at PNs likely to have halos, 
whereas the 388 PNs were collected by random detection. Therefore, the real fraction of 
spherical PNs can be higher than ~ 20 — 30%. The number of PNe known to have halos 
prior to the work of Gorradi et al. (2003) is 25 (R. Gorradi, private communication 2005). 
These amount to ~ 6% of the 388 PNs, rather than the ~ 12% upper limit taken from 
Gorradi et al. (2003). If all PNs are expected to have halos, then the detection fraction of 
non-spherical PNs is ~ 16 times that of spherical PNs. Then the fraction of spherical PNs 
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is ~ 35%. If, on the other hand, only ~ 50 — 60% of PNs have halo (R. Corradi, private 
communication 2005), then the fraction of spherical PNs is again ~ 25%. 

Halos are detected by their surface brightness. In the next section we consider the 
detection fraction from the luminosity of PNs. 


4. ESTIMATING THE UNDETECTED POPULATION BY LUMINOSITY 

SEARCH 

The goal of the present section is to show that it is quite possible that a hidden popula¬ 
tion of spherical PNs exists in detection and surveys based on luminosity rather than surface 
brightness. This will give a lower limit on the hidden spherical population, as spherical PNs 
are even more difficult to detect relative to non-spherical PNs via surface brightness. We will 
build a very simple toy model, then choose typical values for the different physical variables, 
and show that many spherical PNs may have escape detection. This is simply because the 
spherical halo are very dim, hence relative to elliptical and bipolar PNs, spherical PNs are 
much less likely to be detected. We stress that we do not prove the existence of such a 
population. We only show, in line with the discussion in the previous section, that because 
the halos expected to be present around many non-spherical PNs are not detected, many (if 
not most) spherical PNs may escape detection. 

Our goal is to update the statistics of the different morphological classes of PNs and 
their relation to the progenitor type in the binary progenitor theory. Since all large samples 
of PNs which are classified by morphology and to which theory is compared are visible-light 
samples, this study is limited to the discovery of PNs in the visible band. We note however 
that many PNs are bright in the IR, as was already found by IRAS, with typical flux values 
of ~ 1 Jy at 12/im, and ~ 10 — 100 Jy at 60/im (Pottasch et al. 1984). The IR emission 
can be accounted for with models of PN evolution (Volk 1992). PNs are also detected in 
the near IR, e.g., by the 2MASS all-sky infrared survey (Ramos Larios & Phillips 2005). In 
most cases the near-IR emission is dominated by the central star continuum, but extended 
PN are also found (Ramos Larios & Phillips 2005). Although relatively bright in the IR, 
most PNs have been discovered by visible-light emission. Some PNs were discovered by 
their IR emission (e.g.. Beer & Vaughan 1999), but even then [OIII] emission is required 
for dehnite identification. When large enough catalogues of well resolved PNs in the IR 
exists, the same type of studied performed here for visible-light classihcation of PNs can 
be done for IR classihcation. We here note that a hidden population of spherical PNs is 
expected also in samples of PNs discovered by their IR emission, because bipolar PNs, or 
those having symbiotic outhows, are distinguished from other PNs (e.g., M2-9; Ramos Larios 



& Phillips 2005; Schmeja & Kimeswenger 2001). Some bipolar young PNs and proto-PNs 
are particularly bright in the IR, e.g., the Red Rectangle (Miyata 2004), also suggesting that 
a hidden population will exist in future morphological catalogues based on PNs discovered 
by their IR emission. Bright PNs in the near-IR, as another example, were formed by 
progenitors with high mass loss rates (Hora et al. 1999). Such progenitors tend to form 
bipolar PNs. Indeed, most bright PNs in the near-IR are bipolar PNs (e.g.. Latter et al. 
1995, 2000; Hora & Latter 1996; Hora et al. 1999, 2004; Kastner et al. 2001, 2002; Meakin 
et al. 2003). The positive correlation between near-IR emission and bipolar morphology 
(e.g.. Latter et al. 1995) shows that in the near-IR an even larger relative hidden population 
of spherical PNs might exist, adding support to the claim in the present paper. 


4.1. The Toy Model 

Considering the other uncertainties in estimating the progenitor population of PNs, we 
build a toy model to estimate the fraction of the undetected spherical PN population. We 
make the following assumptions. 

1. Spherical PNs are similar in structure to the faint large halo of elliptical and extreme 
elliptical PNs. This assumption is based on the results of Soker (2002), who argues that 
the main difference between spherical and non-spherical PNs with halos, is that the 
progenitors of spherical PNs did not go through the FIW (superwind) phase. Namely, 
the evolution of spherical, elliptical, and extreme-elliptical is similar until the onset 
of a FIW in the final stages of the AGB phase of the progenitors of elliptical and 
extreme-elliptical PNs. This crucial assumption for the model seems to be strongly 
supported by observations, mainly by the similar structure of spherical PNs to that of 
the structure of the halo of elliptical and extreme-elliptical PNs (Soker 2000; 2002). 

2. We assume that for a PN to be discovered, its visible energy flux (observed luminosity) 
should be above some threshold Fq at a considered line, say H(5. Had we used limit 
on the surface brightness rather than on the luminosity, the relative number of un¬ 
detected spherical PNs compared with undetected non-spherical-PNs would be larger, 
strengthening our results (see last paragraph of section 4.3). 

3. We assume that PNs are composed of three parts formed from three mass loss phases. 
A regular AGB wind with a constant mass loss rate of Mh formed the halo. It started 
at time th before present, and ceased at time tg before present. A later FIW with mass 
loss rate of Mg forms the dense shell. It started at time tg before present, and ceased 
at time R before present. Both winds expand with the same speed v. Finally, a faster 
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wind, but of low density, accelerates the inner part of the shell. This forms a dense 
and thin shell, termed rim (Frank et ah 1990). 

4. We assume th ^ i.e., the FIW (superwind) phase is relatively short. 

5. We neglect the time interval between end of the AGB, when we set = 0, and the 
time the PN turns on, namely, when the central star ionizes most of the nebula. This 
is justihed as PNs live for tens of thousands of year, while the post-AGB phase lasts 
few thousands years or less. We also assume that the life time of all PNs is with 
no dependence on the morphological type (no dependence on 77 ; see eq. 1 ). 

6 . We neglect the contribution of the rim to the visible total luminosity, as well as of the 
inner tenuous bubble. We also neglect the change in density structure resulting from 
the propagation of the ionization and shocks fronts (Perinotto et al. 2004). 

7. We assume a spherical PN structure in calculating the luminosity. 

8 . All PNs, with no dependence on morphological type, have the same M^. Following the 
non-detection of halos around many PNs (see previous section), we take this value to 
be as the typical value for long period variables with thin envelopes (Knapp & Morris 
1985), ~ 5 X 10“'^ Mq yr“^. Interestingly, this mass loss rate is Mh ~ O.lL/nc, for an 
AGB stellar luminosity of 3 x lO^L© and wind speed of n = 10 km s“^; c is the speed 
of light. We note that oxygen-rich stars have a median mass loss rate of ~ O.lL/nc 
(Knapp 1986). We take, therefore, the halo mass loss rate to be = 5 x 10 “’^Mq yr“^. 


9. We take the mass loss rates in the two first phases to have a ratio 


Ms 


Mh 


( 1 ) 


Following the discussion in previous sections, we assume that for spherical PNs 77 = 1 
(no FIW in spherical PNs), while 77 > 1 for elliptical and bipolar PNs. 


10. We assume that the mass in the shell that was formed from the FIW, M^, is the 
same for all PNs. As discussed later, taking instead the duration of the FIW, f*, to 
be the same, will make spherical PNs even more difficult to be detected relative to 
non-spherical PNs. 


11. We neglect the different evolution time of central stars having different masses (e.g., 
Blbcker 1995) and take the total time the central star ionizes the nebula to be At* = 
30,000 — 40,000 yr. This assumption, assumption (5), and assumption ( 8 ) above 
are clearly not true for massive progenitors, many of which form bipolar and some 
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which form extreme elliptical PNs. As can be seen from Table (1), this can introduce 
uncertainties of ~ 15—25%, still smaller than some other uncertainties, considering that 
we argue for a hidden population larger by an order of magnitude than the observed 
population. In other words, the parameter most influential in the toy model is rj. 

12. We take most known PNs to be located in a slab, with the galactic plane at its center, 
and we take a constant density of PNs per unit area, apN- Namely, we take a uniform 
two dimensional distribution of PNs, neglecting the high density toward the galactic 
plane. We also take most PNs to be in the distance range of 1 kpc ~ Dmin < D < 
Dmax ~ 5kpc. 


Despite being strong, these assumptions capture the basic factors involved in a possible 
undetected population of spherical PNs. 

The luminosity of the nebula L is the sum of the halo and shell luminosity. The halo 
is bounded between = vtg and Vh = vth, and its density is ph = M/j/(47rr^n); a similar 
expression holds for the shell bounded by = vtr and r* = vtg- The total mass in the halo 
is Mh = Mh{th — ts), and that in the shell Mg = Mgitg — A). We find 


MMi-D+Mdl-l 

n rgj \rs Th 


MgMg ^ Ml{th - ts) 






thts 


( 2 ) 


where dV = Aiir'^dr is a volume element, and A; is a constant such that I = kp^ is the 
emissivity (power per unit volume). At a temperature of 10^ K we find for Hj3 that k = 
2.6 X 10^^ erg cm^ s“^ g~^. Using our assumption 3> U, and the relation tg = tr + Mg/Mg, 
equation (2) becomes 


L 


k 

47j-p3 




Using the definition of p (eq. 1), defining 


Mh 


X 10- 




M, 


-1 


0.2Mo / \ 5 X lO-W© yr 


.-1 


yr, 


( 3 ) 


( 4 ) 


^ k Ml^ 

Ativ'^ t ’ 


and 


L 


( 5 ) 
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we can write equation (3) as 


L — Lg 1 + 77 


r 

'Tr 


1 + 


1 r 
rjtr 


-1 


r 


( 6 ) 


By our assumption (2), for detecting a planetary nebulae F = L/{A7rD'^) > Fq, where 
D is the distance to the PN. We write this condition as 


F - FgD,^ ( 1 + ?7— ( 1 + — > Fo, 


kpc 


1 r 


rjtr 


-1 


( 7 ) 


where -Dkpc is the distance to the PN in units of kpc, and 


Fg{H(3) = 


47r(l kpc) 


~ 10 


-15 


Mu 


Mg 


-1 


5 X lO-Wo yr-i / \ 0.2M£ 


© 


X 


-3 


erg s ^ cm ^ 


-10 km s“^ j o - - (8) 

Solving for tr, we hnd the time during which the PN can be detected. Equation (7) can be 
cast into a quadratic equation 






( 9 ) 


Neglecting the post-AGB evolution time (assumption 5 above), the time span during 
which a PN at distance D with enhanced mass loss rate factor rj can be detected, is given by 
the solution of the last equation. We scale the time with r, and marked it by At{ri, D) = tr/r. 
However, there is another condition on At. In the treatment above we have neglected the 
fading of the ionizing central star. The condition we should add is that the PN life time will 
not exceed the time it takes the central star to fade, which itself depends on the central star 
mass. Ignoring this in our toy model, assumption (11) above, we add therefore the condition 

At{T], D) <Atc = — ^ 0.1, (10) 

T 

where the values used come from equation (4) and assumption (11). The last condition 
implies that for a distance D < Ddr}) we should set At{ri,D) = Ate from equation (10), 
rather than taking the solution of equation (9). Solving equation (9) and using condition 
(10), gives the time span during which the PN can be detected 


At{r], D) = min 


1 ( 
2\Fo 



-2 




1/2 



( 11 ) 
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The case 77 = 1 has a simple solution 


At{ri = l,D)= min \ , At, 

-To 


( 12 ) 


The critical distance Dc below which Ate should be used is found by substituting tr/r = 
Ate in equation (9). This gives 


De = 


Fs Ate + T] 
Fn Ate 


Ate + 


V 


-1 


1/2 


(1 + lOf]) 1 + 


10 

V 


-1 


1/2 


(13) 


where in the second equality we substituted the typical numerical values used here, Fs/F q = 
0.1 (see next subsection), and Ate = 0 . 1 . 


4.2. Numerical Values 

To estimate the value of rj we start with the density prohle of the PN NGC 6826 as 
given by Plait & Sober (1990). This PN has a large spherical Halo, with an elliptical inner 
structure, similar to several other PNs (Corradi et ah 2003). For a constant mass loss rate 
the density decreases with radius as R~^. In the boundary between the inner elliptical region 
and the outer spherical halo the density gradient is much steeper. In going from ~ 12 arcsec 
to ~ 18 arcsec the density decreases by a factor of ~ 5. Comparing to (18/12)^ = 2.25, 
we hnd that the increase in mass loss rate was r] ~ 2.5. In a study of the proto-PN HD 
56126 Meixner et ah (2004) hnd that the mass loss rate during the superwind (FIW), which 
lasted 840 years, was 77 = 6 times higher than that in the previous mass loss episode, 
M = 3 X 1 O“^M 0 yr“^ versus M = 5.1 x IO^^Mq yr“F We therefore take 77 = 4 to be our 
standard value. 

The value of the detection hux limit is very hard to give, because PNs are discovered 
by many different telescopes and surveys. Our approach is to examine the histogram of H(3 
hux of PNs from the catalogue of Acker et ah (1992; their hg. 9a). The lower detection hux 
is ~ 10“^^ erg s“^ cm“^, with the majority of PNs with F > 10“^"^ erg s“^ cm“^, which we 
take to be the detection limit Fq. Therefore, our toy-model value is Fs/Fq = 0.1, where Fg 
is by equation ( 8 ). 


4.3. Detection Probability 

The chance to detect a PN of a morphological type which we characterize by the en¬ 
hanced mass loss rate factor 77 , at a distance D, is the fraction of time for which it obeys 
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equation (7), which solution is given by equation (11), 


P^d(-D) oc 


tdv) 


N,{D), 


(14) 


where tii'f]) is the life time of the PN, assumed not to depend on r] in our toy model, 
and Nj^{D) = 27iap-^DdD is the number of PNs between D and D + dD, and we used 
our assumption (12) of constant PN density per unit area in the disk. Normalizing the 
probability under our assumption of PNs slab distribution, gives 


»Du 


Pv = 


7r(^Lx - ^mm)c^PN 


27rZl(TpN 


'D„ 


tdv) 


dD 


(15) 


As we are interested only in relative probabilities among the different morphological types, 
marked by rj in our toy model, we can cancel common factors, such as the life time II- The 
relevant probability is then 

P 1 

= p(^ 11) = p®) 


where with the aid of equation (12), we hnd 


cDj, 


Pi = / DAtiv = 1, D)dD = ^ In ^ 

n 0 ^min 


~ 0.1 InPp 


'Dn 


( 17 ) 


where in the second equality we substitute our standard values of Pmin = 1, and P^/Pq = 0.1. 
For 7] ^ 1 such that Dc > Pmax, 

Pr,{D, > Pmax) = ^ DAtJD = ^ ^ 

where in the second equality we substitute our standard values of Ate = 0.1, Pmin = 1, 
Pmax = 5, and P,/Po = 0.1. For P^ax = 4 and 6, we obtain (for Dc > Pmax) Pr», = 5.4, 
and 9.8, respectively. In Figure 1 we draw Pr^ as function of r] for the above parameters. 
A good £t to the line for the above parameters is Puriiv < 6.4) = 0.00682?7^ — 0.10864?7^ + 
0.360827]^ + 1.67047] — 0.95117, while for rj > 6.4, as we found above, Pr,,(7] > 6.4) = 7.46. 

When examining the results presented in Figure 1, one should recall the purpose of 
this section. We here do not prove and do not show that there is a large population of 
undetected spherical PNs. This claim (see previous section) is based on the small detection 
fraction of spherical halos around elliptical and bipolar PNs (Corradi et ah 2003), and our 
basic postulate that most spherical PNs are similar in structure to spherical halos around 
non-spherical PNs. In the present section we limit ourselves to show that it is indeed possible 
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that most spherical PNs are below detection limits of previous observations. In other words, 
in this section we show that there is no contradiction between the claim that many spherical 
PNs are not detected, and the known properties of PNs and halos of PNs, under our basic 
postulate. For that, in this section we built a toy model, and took reasonable values for 
the physical parameters. Figure 1 shows just that, and no more. Namely, for our favored 
toy model value for the mass loss enhancement factor r] = A, the detection probability of 
spherical PNs is an order of magnitude below that of non-spherical PNs. This value is only 
suggestive, because of the nature of the toy model, but it compatible with the spherical-PN 
detection fraction of ~ 10% found in section 3. 

The toy model actually overestimates the detection probability of spherical PNs relative 
to that of non-spherical PNs. Our assumption (10) that all PNs have the same shell mass 
Ms, was implemented by scaling the time with the same r (eq. 4). Lower FIW mass loss rate 
(lower f]), implies lower density, but wider shells. As the luminosity is proportional to density 
square times the volume, the luminosity of the shell is proportional to rj. Had we taken the 
duration of the FIW that formed the shell to be the same for all PNs, which implies the 
same shell’s volume for all PNs, the shell’s luminosity would have been proportional to r]'^. 
Namely, the detection probability would have increases faster for higher mass loss rate. The 
same argument goes for a calculation to hud the detection probability by surface brightness. 
Since the shell has the highest surface brightness, the detection of PNs is by the shell, 
with negligible contribution from the halo (the relative contribution of a large halo to the 
luminosity is more signihcant than to the surface brightness). The surface brightness of the 
shell goes as ~ ? 7 ^, and not as rj. For these considerations we have build our toy model to be in 
the safe side, in the sense that using other reasonable assumption (like considering detection 
by surface brightness) will only strengthen our claim for hidden population of spherical PNs. 


5. SUMMARY 

The goal of the present paper is to update the statistical distribution of different mor¬ 
phological types of PNs and their progenitors, in the frame of the binary-shaping model. 
This update is motivated by, and uses, recent results on the distribution in mass and orbital 
separation of extra-solar planets, and the new work of Corradi et al. (2003) on halos around 
PNs. We make two assumptions: 

1. We take the common view that the axisymmetric (including point symmetric) structure 
of PNs result from binary interaction; the companion can be stellar or substellar. 

2. We assume that spherical PNs (according to the dehnition of Soker 1997 which is used 
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here) are similar in structure to the spherical halos around non-spherical PNs. This 
assumption is based on the results of Soker (2002). 

Only the spherical morphological type is studied here. We conclude that most probably there 
is a large population of spherical PNs, much more than non-spherical PNs, that are too faint 
to be detected. This leads to a new estimate regarding the distribution of PN progenitors 
and different PN morphology types (summarized in sections 1 and 2). The new estimates 
compared with previous ones are summarized in Table 1. 

To reach the conclusion that there is a hidden spherical PN population, we hrst examine 
Corradi’s et al. (2003) observations and analysis of halos around PNs. Whereas all PNs are 
expected to have a faint halo, the detection fraction is much lower; taking the PNs listed 
by Corradi et al. (2003) to have detectable halos, and compare with all PNs, we hnd that 
only ~ 12% of all PNs have detectable halos. Using our second assumption above, this 
suggest that the detection fraction of non-spherical PNs is ~ 8 times higher than that of 
spherical PNs. Because many of the PNs listed by Corradi at al. (2003) were a target 
of specihc (deeper) observation, the average (among all large surveys) detection fraction of 
non-spherical PNs can be ~ 10 — 15 times higher than that of spherical PNs. We therefore 
claim that while ~ 3% of all well resolved PNs are spherical (Soker 1997), the true fraction 
of spherical PNs among all PNs is likely to be ~ 30%. The recent detection of a spherical 
PN at a distance of only ~ 0.5 kpc (Frew & Parker 2006) show that spherical and faint 
PN can indeed escape detection. This conclusion regarding a hidden population is not to be 
confused with the claim of Manchado et al. (2000) for ~ 25% observed spherical PNs; they 
simply use a different dehnition (which we don’t accept; see section 1) for spherical PNs. 

We then (section 4) built a toy model to show that such a large (relative to that of 
non-spherical PNs) hidden population of spherical PNs might be compatible with present 
sensitivity of observations. We stress here that in section 4 we do not determine or estimate 
the detection fraction of spherical PNs. The uncertainties in many physical parameters (in 
particular in the progenitor’s mass loss history) are too large to allow such an estimate. 
However, the toy model shows that the claim for a large hidden spherical PN population 
(section 3) is compatible with what we know about PNs. 

One interesting conclusion from this paper concerns the fraction of PNs that were shaped 
by planetary systems. As seen from Table 1, the required fraction of PNs that were shaped 
by planets was reduced substantially from that estimated by Soker (1997). The stellar 
population synthesis studies yield a relatively secure fraction for the PN shaped by stellar 
companions. The three stellar evolutionary routes add up to < 60%. The rest of PN progen¬ 
itors either yield a spherical PNs or were interacting with substellar objects to form elliptical 
PNs. In the past, most of these progenitors were assumed to have substellar companions that 
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influenced the mass loss geometry. The present study reduces this number. This is in accord 
with the hnding that the fraction of Sun-like stars that have close and massive planets or 
brown dwarfs around them is < 25% (see section 2). We actually argue that the number of 
AGB stars that formed spherical PNs is larger than those having substellar companions that 
influenced their mass loss process. Our best estimate in the present study is that the hidden 
population of spherical PNs, although large, is too small to account for all PNs not shaped 
by stellar companions; hence shaping by planetary system seemed to be required in some 
cases. However, we cannot completely rule out the possibility that the hidden population of 
spherical PNs is larger, and that practically no PN was shaped by planets (some may have 
been shaped by brown dwarfs). 

The hidden population of spherical PNs, if exist, will not change much the comparison 
between the total number of PNs, and the number expected from stellar counts. This is 
because the proposed hidden population adds only ~ 30% to the total number of PNs, a 
number much smaller than the uncertainties in both population synthesis and in estimating 
the total number of PNs in the galaxy. This is the place to stress again that our finding 
regarding the hidden population of spherical should be considered in the context of the 
general estimate of PN progenitor population presented in Table 1. For example, the hidden 
population is compatible with the hnding that the fraction of stars having brown dwarfs or 
massive enough planets to shape the descendant PNs is ~ 25%, but only a fraction of them 
(estimated here to be ~ 15% will form PNs (see section 2). 

The main prediction of this paper is that very deep observations will detect more spheri¬ 
cal PNs. This is not easy. After all, instead of the 14 spherical PNs classihed by Soker (1997), 
we expect that for the same sample of PNs there will be ~ 150 spherical PNs, which is a 
small number for deep observations of particular patches of the sky (they will not be highly 
concentrated in the galactic plane). In addition, they are expected to be very faint, hence 
to be detected they require long observations. A survey such as the AAO/UKST Ha one 
(Parker et ah 2003), might discover higher percentage of spherical PNs. Indeed, a circular 
shell around an evolved star that was detected by this survey (Cohen et al. 2005). However, 
this survey is based on surface brightness. As discussed at the end of the previous section, 
in surface brightness surveys spherical PNs are even less likely to be discovered compared 
with non-spherical PNs. Therefore, it is not clear whether the proposed hidden spherical 
PNs population will be revealed by that search. 


We thank Letizia Stanghellini, Romano Corradi, and Joel Kastner for useful comments. 
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Table 1. 


BINARY INTERACTION AND MORPHOLOGY OF DESCENDANT PNs 


Binary interaction 

Major PN type 

1997 

Percentage 

2001 

This paper 

Close stellar companion 
outside envelope 

Bipolar 

~ 11 

~ 15 

~ 15 

Stellar companion in 
a common envelope 

Extreme 

Elliptical 

~ 23 

~ 25 

~ 25 

Substellar companion 
in a common envelope 

Elliptical 

~ 56 

~ 35 

~ 15 

Wide stellar companion 

Elliptical+jets 

— 

~ 15 

~ 15 

No interaction 

Spherical 

~ 10 

~ 10 

~ 30 


Note. — Percentage of the different morphological types and their progenitors 
in the binary shaping model. The second column gives the major, but not the 
sole, PN descendant from the mode of binary interaction. Percentage are of all 
PNs, and not of the observed PNs only. As discussed in the paper, most spherical 
PN are too faint to be detected. The percentage are from Sober (1997), Sober 
(2001a,b), and present paper. 
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Fig. 1.— The probability of detecting a PN in the distance range 1 — 5 kpc, as a function of 
the enhanced mass loss rate factor rj (eq. 16). In the toy model, rj is the ratio between mass 
loss rate during the hnal intensive wind (superwind) on the AGB and the regular AGB wind 
(eq. 1). The probability is given relative to the probability of detecting a PN with no 
enhanced mass loss rate, which is assumed to be spherical. 




